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We consider asymptotics of the velocity derivatives probability distribution functions (PDFs) in 
Burgers' turbulence. We argue that in the forced case the same power laws as in the decaying case 
are realized for an infinite system. 



PACS numbers: 47.27.Gs, 02.50.Ey, 05.60.-k 



The Burgers equation describes an evolution of weakly 
non-linear one-dimensional acoustic waves in the refer- 
ence frame moving with the sound velocity ffl . The equa- 
tion is written for the velocity u 



dtu + ud x u — va^u = f . 



(1) 



where v is the viscosity coefficient and / represents an 
external influence pumping the waves. The most inter- 
esting feature of this equation is that it describes shock 
formation. 

Recently statistical properties of solutions of the Burg- 
ers equation ([!]) have attracted a considerable attention, 
since it is a simple model which could demonstrate some 
features of developed turbulence [0 . Then one assumes a 
large value of the Reynolds number Re = VL/v, where V 
and L are characteristic velocity and length of the large- 
scale motion. One can treat decaying (/ = 0) or forced 
turbulence. In the decaying case statistics of u is time- 
dependent, whereas it is expected to be stationary, if the 
pumping force / is a process statistically homogeneous 
in space and time. We consider both possibilities. 

An object, mainly treated in papers devoted to Burgers 
turbulence, is the PDF "Pi(£) of the velocity derivative 



This PDF is strongly asymmetric. The right tail 



of the PDF decays fast in accordance with the inviscid 
law InT^i oc — £ 3 The inviscid left tail entailed a 

controversy. Different values of the exponent a for this 
tail V\ tx |£|~ Q were proposed in the works |4|Jt^[ll|l, see 
also the discussion in [ fn^|l3| . Below we confirm the va- 
lidity of the answer a = 7/2, claimed in |) 11 1, for the 
forced case in an infinite system. 

We assume that the system size C is much larger than 
L. Then the PDF of the velocity derivative £„ = 
can be written as the space average 



V n {t^ n )=C- 1 / dx5[£ n 



- „(») 



(t,x)} 



(2) 



The expression is the starting point of our consideration. 
It is argued in the works P Jll| , |l4j that the inviscid left 
tail of V\{£) is related to the pre-shock regions. These 
vicinities possess some universal properties enabling to 
establish the PDF (@Y 



We begin with the decaying case when / = in Eq. 
(0), assuming a smooth initial condition u (x) at t = 0. 
This case was examined by Bee and Frisch pJJ. Below, 
we reproduce their arguments in a slightly different lan- 
guage. 

It is well known that Lagrangian trajectories coming 
to the shock creation points start at t = from points 



with zero second derivati 



= 0. First terms of the uq 



expansion near such point X are 

m ~ v + a(x - X) + -(x - X) 3 . 

6 



(3) 



The shock is produced provided a < and s > 0. Clearly, 
the contribution to the integral (||) related to the shock 
is independent of X and v. Therefore the integral (^) is 
rewritten as 



P n (t,S n ) = — J dads TV, s) J dx5[& 



,(»)! 



(4) 



where V(a,s) is the PDF of the parameters a,s in the 
points Xi and Af is their number. The integral over x in 
Eq. (|J) is taken over a vicinity of the shock formation 
point related to the initial velocity profile (||) . To be pre- 
cise, a number of Lagrangian trajectories, starting from 
the points Xi, are "eaten" by shocks, already formed, be- 
fore coming to the shock creation point. Such points Xi 
have to be excluded from the consideration, that leads to 
a ^-dependence of both Af and V(o, s) in Eq. (|J). 

The shock is formed at t = | cr | 1 . If \a\t < 1 then we 
can use the characteristics method, which enables one to 
write u(t, x) in the parametric form 



x = y + ut , 



One can explicitly take the integral over x in Eq. (Q) 
passing to integration over y. For large £ corresponding 
to 1 — \a\t <C 1 we get 



dx6[£ 



|£|3 t 5/2 



1 



-1/2 



with the restriction 1 — \cr\t < (t|£|) 1 (otherwise the 
integral is zero). Then the integration over a leads to 



1 



^-Z^/^-'-V). (5) 

Thus we reproduce the answer obtained in the 
works [0Jll||- Analogously, the inviscid Khanin law 
|^|-(3n^i)7l3n-i) ( wr itten in Q) can be obtained. 

Now we proceed to investigating the forced Burgers 
turbulence. We consider the force / continuously dis- 
tributed over time and space and having a homogeneous 
statistics. Then the statistics of £ n is also homogeneous 
in time. We assume that / is smooth both in space and 
time. 

As previously, we are interested in vicinities of the 
shock formation points. Let us consider a prehistory 
of the u evolution during the time t before the shock 
creation point, where L/V ^> t ^> The first in- 

equality enables one to neglect the pumping force during 
the evolution. The second condition ensures the validity 
of the expansion (|j|) for examining the shock formation. 
Thus, we return to the expression (||), where V(a,s) is 
the PDF of a = u' and s = u'" in points where u" = 0, 
u' < 0, u'" > 0, this PDF is stationary in the pumping 
case. However, there is a problem related to an explicit 
dependence of the answer on the auxiliary time t. Let us 
demonstrate, that such dependence is really absent. 

The inequalities imposed on the evolution time t mean 
that the values 1 <C crL/V, sL 3 /V are relevant in the 
integral (||). We assume also crL/V, sL 3 /V < Re 1/2 . In 
this range one can neglect both pumping and the viscous 
term in Eq. ([[]). Then we get from Eq. (Q) a closed sys- 
tem of dynamical equations for a = v! and s = vl" (at 
the condition u" = 0), leading to the following equation 
for their PDF 

d t V = d a (a 2 V) + d s (AasV) , (6) 

For the stationary case the time derivative in Eq. (^) is 
zero and we get 

V = a- 2 s- 1 T(s/<7 i ). (7) 

Of course, the function F{z) in M) knows about pump- 
ing. Substituting the expression (^) into Eq. (^) we get 



Vl{0 ~TW*I % Hz) - (8) 

We observe that the auxiliary time t drops from the 
expression. Thus, we reproduce the law |£|~ 7 / 2 for 
the pumping case. The same procedure gives the law 
|£j-(3n+4)/(3n-i) f or high derivatives. 
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